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Abstract. The Cauchy problem for the three-dimensional compressible flow of nematic 
liquid crystals is considered. Existence and uniqueness of the global strong solution 
are established in critical Besov spaces provided that the initial datum is close to an 
equilibrium state (1,0, d) with a constant vector d £ S 2 . The global existence result is 
proved via the local well-posedness and uniform estimates for proper linearized systems 
with convective terms. 

1. Introduction 

Liquid crystals are substances that exhibit a phase of matter that has properties between 
those of a conventional liquid, and those of a solid crystal. The three-dimensional flow 
of nematic liquid crystals can be governed by the following system of partial differential 
equations [U [IT] : 

|? + div(pu) = (1.1a) 
d^"- 1 + div(pu (g) u) - /iAu -(fi+ A)Vdivu + VP(p) (1.1b) 
= -£div (vdQ Vd- i|Vd| 2 J^ (1.1c) 

— + u-Vd = 6?(Ad + |Vd| 2 d), (l.ld) 

where p G M is the density function of the fluid, u G R 3 is the velocity and d £ S 2 represents 
the director field for the averaged macroscopic molecular orientations. The scalar function 
P 6 R is the pressure, which is an increasing and convex function in p. They all depend 
on the spatial variable x = (x±, X2, X3) E M 3 and the time variable t > 0. The constants p, 
and A are shear viscosity and the bulk viscosity coefficients of the fluid respectively that 
satisfy the physical assumptions p > 0, 2p + 3X > 0. £ > 0, > stand for the competition 
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between the kinetic energy and the potential energy, and the microscopic elastic relaxation 
time (or the Debroah number) for the molecular orientation field, respectively. The symbol 
(g> denotes the Kronecker tensor product such that utgiu = (ujUj)i<ij<3 and the Vd©Vd 
denotes a matrix whose ij-th entry (1 < i,j < 3) is d Xi d ■ d x d. I is the 3x3 identity 
matrix. To complete the system the initial data are given by 

p\t=o = po(x), u\ t =o = u (x), d| t =o = d (x), with d £ S 2 . (1.2) 

Roughly speaking, the system (jl.ip is a coupling between the compressible Navier-Stokes 
equations and a transported heat flow of harmonic maps into S 2 . It is a macroscopic 
continuum description of the evolution for the liquid crystals of nematic type under the 
influence of both the flow field u, and the macroscopic description of the microscopic 
orientation configurations d of (rod-like) liquid crystals. 

The hydrodynamic theory of liquid crystals in nematic case has been established by Er- 
icksen and Leslie, see [71151 1151 116|. Since then, the mathematical theory is still progressing 
and the study of the full Ericksen-Leslie model presents relevant mathematical difficul- 
ties. When the density is a constant, the system (jl.ip becomes an incompressible model. 
In [T7], Lin introduced a simplification of the general Ericksen-Leslie system that keeps 
many of the mathematical difficulties of the original system by using a Ginzburg-Landau 
approximation to relax the nonlinear constraint d £ S 2 . Later in [T5] , Lin and Liu showed 
the global existence of weak solutions and smooth solutions for that approximation system. 
For more results on the approximation system, we refer to |191 122"1 124] , Recently, Hong 
[llj and Lin-Lin-Wang [20] showed independently the global existence of weak solution 
of an incompressible model of system (jl.ip in two dimensional space. Moreover, in [20] . 
the regularity of solutions except for a countable set of singularities whose projection on 
the time axis is a finite set had been obtained. In the recent work [23] . Wang established 
a global well-posedness theory for the incompressible liquid crystals for rough initial data, 
provided that 

lluollsA/O- 1 + [do\BMO < £o 

for some eq > 0. Note that the relationship between BMO~ l and is (see [2]) 

m ^ L 3 ^ Bp£* BMO' 1 , for3<p<oo. 

Concerning the compressible case, local existence of unique strong solutions of (jl.ip was 
proved provided that the initial data pq, Uq, dg are sufficiently regular and satisfy a natural 
compatibility condition in a recent work [13] . A criterion for possible breakdown of such 
a local strong solution at finite time was given in terms of blow up of the L°°-norms of p 
and Vd. In |12j , an alternative blow-up criteria was derived in terms of the L°°-norms of 
Vu and Vd. The global existence of weak solutions to (jl.ip with large initial data is still 
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an outstanding open problem for high dimensions. By so far, only results in one space 
dimension have been obtained, for instance, we refer to [5j[6]. 

In this paper, we are interested in the existence and uniqueness of global strong solutions 
to the Cauchy problem of (jl.ip with initial data (jl.2p in the three dimensional space. It is 
difficult to find a functional space such that the system (jl.ip - (jl.2p is well-posed globally 
in time. To this end, we notice that the system (jl.ip is invariant under the following 
transformations 

p = p(Z 2 i,Zx), u = lu(l 2 t,lx), d = d(l 2 t,lx) (1.3) 

with the modification of the pressure P = l 2 P. A critical space is a space in which the 
norm is invariant under the scaling 

(e,f,g)(z) = (e(lx),lf(lx),g(lx)). 

The scaling invariance (jl.3p reminds us of a similar property of three dimensional in- 
compressible Navier-Stokes equations, which provides a well-known global existence of 
solutions with small data in the homogeneous Sobolev space H 2 (see [H]). Motivated by 
this observation, we aim at a global well-posedness of the system (|l.ip - ([1.2|) with small 
initial data in a functional framework where the function space for the velocity u is similar 
to H*. According to the scaling (II. 3D . the regularity of the density p and the director 
field d is one order higher than that of the velocity u, and hence a function space which 

■ 3 ■ 3 

is similar to would be a candidate. Unfortunately, the function space does not 

■ 3 

turn to be a good candidate for p and d, since bounds of the director field d in H 2 do not 
automatically imply the L°° bound of d. To overcome this difficulty, inspired by [3] for the 

compressible Navier-Stokes equations, it seems more natural to work in the framework 

3 3 

of homogeneous Besov space B 21 srnce ^21 * s continuously embedded into L°°. Further- 
more, as in [3], the different dissipative mechanisms of low frequencies and high frequencies 

_ 1 3 3 1 
inspire us to deal with p and d in B 2 2 1 2 = B% 1 n B 2 1 ■ 

The rest of this paper is organized as follows. In Section 2, we will give a series of 
fundamental properties of the Besov's spaces. In Section 3, we reformulate the system 
(|l.ip - (|1.2j) and state our main result. The main goal of Section 4 is to prove uniform 
estimates for linearized systems to (jl.ip . while the global existence is obtained in Section 
5. In Section 6, the uniqueness of global strong solution is verified. 

2. Preliminaries 

Throughout this paper, we use C for a generic constant, and denote A < CB by A < B. 
The notation A ~ B means that A < B and B < A. Also we use (a q ) q ^z to denote a 
sequence such that X^ez ^ < 1- (f\d) denotes the inner product of two functions f,g in 
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L 2 (R N ). The standard summation notation over the repeated index will be adopted in 
the remaining part of this paper. 

In order to state our existence result, we introduce some functional spaces and explain 
the notations. Let T > 0, r E [0, oo] and X be a Banach space. We denote by Ai(0, T; X) 
the set of measurable functions on (0, T) valued in X. For / E M(0,T; X), we define 

i 



\l°°(X) = supess re(0iT) ||/(r)||x. 

Denote L r (0,T;X) = {/ € M(0, T; X)\ ||/|| £ r {x) < oo}. If T = oo, we denote by 
L r (R + ;X) and ||/||.z/-pn the corresponding spaces and norms, respectively. C([0,T], X) 
(or C(R + ,X)) stands for the set of continuous X- valued functions on [0, T] (resp. M + ) 
while Cfe(]R + ;X) is the set of bounded continuous X-valued functions. For a G (0,1), 
C a ([0,T]; X) (or C°(M + ;X)) stands for the set of Holder continuous functions in time 
with order a, i.e., for every i, s in [0,T] (resp. we have 

\\f(t)-f(s)\\x<\t-sr. 

As in 0, we introduce a function ip E C°°(R N ), supported in C = {£ E R N : § < |£| < 
t^} and such that 



Let J 7 be the Fourier transform. Denoting J 7 by h, we define the dyadic blocks as 
follows: 



p<q-l 



A q f = V(2"*Z>)/ = 2^ / /i(2«y)/(x - and S q f = V A p /. 

Then the formal decomposition 

is called homogeneous Littlewood-Paley decomposition. 
For s E K and / E S'(R N ), we denote 



de_f 
B s ~ 



.qez 



When p = 2 and r = 1, we denote || • \\p, s by I ■ ||s s - The definition of the homoge- 

p,r 

neous Besov space is built on the homogeneous Littlewood-Paley decomposition (cf. j3j 
Definition 1.2]) 
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Definition 2.1. Let s£l, and m = — [y + 1 — s] . If m < 0, we set 

= | / G ^(M^III/Hb. < oo and / = £ A,/ in S'(l^) I . 

If m > 0, we denote by V m the set of N variables polynomials of degree < m and define 

B°=\f€ S'(R N )/r m \\\f\\ B s < oo and / = ^A q f in S'(R N )/V n 

Functions in the homogeneous Besov space B s has many good properties (see [U Propo- 
sition 2.5]): 

Proposition 2.1. The following properties hold: 

• Density: the set Cq° is dense in B s if \s\ < 

• Derivation: \\f\\B s ~ ||V/||b«-i; 

• Fractional derivation: let A = \J — A and er G M; f/ien i/ie operator A a is an 
isomorphism from B s to B s ~ a ; 

• Algebraic properties: for s > 0, B s n L°° is an algebra; 

• Interpolation: (B Bl , B"*)o tl = B 6si+ ^~ 6 > 2 , for Sl ,s 2 G M and G (0, 1). 

To deal with functions with different regularities for high frequencies and low frequen- 
cies, motivated by [3], it is more effective to work in hybrid Besov spaces. We remark 
that using hybrid Besov spaces has been crucial for proving global well-posedness for 
compressible Navier-Stokes equations in critical spaces (see [HE]). 

Definition 2.2. Let s,t£R. We set 

||/||^, t =^2^||A g /|| L2 +^2^||A g /|| L2 . 

q<0 q>0 

For m = — + 1 — s] , we define 

B^ = {feS'(R N )\\\f\\ §s , t <oo} ifm<0, 
B 8 * = {f e S\R N )/r m \\\f\\Bs, t < oo} if m>0, 

Remark 2.1. Some remarks about the hybrid Besov spaces are in order: 

• B s ' s = B s ; 

• Ifs < t, thenf? 5 '* = B S C\BK Otherwise, B s ' 1 = B s +B l . In particular, £? s 'f > L°° 
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The space B°' s coincides with the usual nonhomogeneous Besov space 
/ G S\R N )\\\ X (D)f\\ L2 + ^2^||AJ|| L2 < oo 

q>0 
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where x (0 = 1 - E q >o H^Ol 

• If si < s 2 and t x > t 2) then 5 s1 '* 1 ^ B 3 ^ . 

We have the following properties for the product in hybrid Besov spaces (see pQ): 
Proposition 2.2. • For all s\,S2 > 0, 

1 \fg 1 1 =2 ^ II/IUHMIbsi^ + ll/llfl«i,»a ; 

• For all si, s 2 < y suc/i i/iaf minjsi + ti, S2 + £2} > 0, 

||/3||^ S1+S2 _^ ,t 1+ t 2 -f < ll/llsn.^llS'llBti.ta- 

Throughout this paper, the following estimates for the convection terms arising in the 
linearized systems will be used frequently (cf. Lemma 5.1]). 

Lemma 2.1. Let G be an homogeneous smooth function of degree m. Suppose — ^ < 
Si,t < 1 + y for i = 1,2. Then the following three inequalities hold true: 

\(G(D)A q (u-Vf)\G(D)A q f)\ 

< Ca q 2-^ n ^)-m )||u|| sl+ ^||/y si , S2 \\G{D)\f\\ L2 , (2,1) 
\{G{D)A q (n ■ Vf)\A q g) + (A,(u • V ff )|G(D)A 3 /)| 

< CoalluH Rl+4 (2-* HG^A,/!!^ + 2 "«C^(«)-m) ||/|| ^ 1| A g5 || L2 ) , (2 ' 2) 



where ^2 q€ %oi q < 1 and C is a universal constant that only depends on Si,t,N. The 
notation (j) s,t (q) means that for q £ Z, 



s, i/ g < 0, 
i, i/ g >l. 



Remark 2.2. The above lemma looks slightly different from the original one (3] Lemma 
5.1]. Indeed, (2.2) and (2.3) are corresponding to (41) and (42) in [3], respectively, up to 
a change with the notation </> s '*. 



3. Reformulation of the Original System (|l.l|) and Main Result 

In this section, we first reformulate the original system (jl.ip into a different form and 
then we state our main result on the global existence of strong solutions. We simply set 
6 = ^ = 1 since their sizes do not play any role in our analysis. For s £ R, we denote 

Using the idea in [3], we decompose the velocity field into a compressible part and an 
incompressible part. Let 

h = A _1 divu 
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and 

= A _1 curiu, with (curlu)^ = d x u l — d Xi u 3 . 
Owing to the identity A = Vdiv — curlcurl, we have the decomposition 

u = -A _1 V/i + A -1 curin, 

which implies that u can be recovered from the information of h and Q. Denote 



de f 

A = fiA + (A + p)Vdiv, 



and 



(3.2) 



N = -u • Vu - t — l^u - idiv ( Vd Vd - -|Vd| 2 / ) . (3.1) 

p p V 2 / 

Applying A _1 div and A -1 curl to the moment equation in (jl.ip respectively, we obtain 
that 

' d t h - vAh = A-Miv (j\f - ±VP(p) N 

d t Q - /iAfi = A-^urW, 

where u = 2ji + A. Since fi > and 2/i + 3A > 0, we have ^ > |/i > 0. In the second 
equation of (|3.2p . we have used the fact that 

curl (^-VP(p^j = -curlVP(p) + curl ^-J x VP(p) = 0. 

The advantage of the above reformulation is to get rid of the pressure for f2, the incom- 
pressible part of the velocity, while we still keep all information of the velocity field u. 
For the simplicity of our presentation, our proof focuses on the case: 

P(p) = l p 2 . The 

general barotropic case (P(p) is an increasing convex function of p) can be verified by a 
slight modification of the argument below. 

In this paper, we shall prove the existence of global strong solution for initial datum 
that is close to an equilibrium state (1,0, d) with a constant vector d £ S 2 . The result is 
valid for any positive constant density p and we take p = 1 just for simplicity. Keeping 
(|3.2p in mind, it is convenient to reformulate the original system (jl.ip into a new system 
in terms of p, h, Q and d 

d t (p - 1) + u • V(p - 1) + Ah = -(p - l)divu, (3.3a) 

dth - uAh - A{p - 1) = A _1 div7V, (3.3b) 

d t tt - pAn = A^curLV, (3.3c) 
3d 

— +u-Vd- Ad = |Vd| 2 d, (3.3d) 
ot 

subject to initial conditions 

p\t=o = Po(x), h\ t=0 = h (x) = A _1 divu (x), 
£l\t=o = = A~ 1 curlu (x), d|t =0 = d (x). 



(3.4) 
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Let us now introduce the functional space that appears in the global existence theorem. 

Definition 3.1. For T > 0, and s G R, we denote 

58£ = {(e,f,g) G (l 1 (0,T;B s + 1 ' 8 )nC([0 ) T];B s - 1 ' s )) 
x (L 1 (0,T;S s+1 )nC([0,T];B s - 1 )) 3 
x (^(O, T; j5 s+1 > s + 2 ) n C([0, T]; S 5 - 1 ' 5 )) 3 } 

and 

ll( e ,f ; g)ll<Bf, = ll e IL°?(B«-i>*) + W^W L^iB"' 1 ) + l|glL°?(B<-M) 

+ ll e llz4(Bs+M) + H f llz4(Bs+i) + ||g|| L i,(B S +i, s +2)- 

We use the notation 33 s if T = +oo by changing the interval [0, T] into [0,oo) in the 
definition above. 

The main result of this paper is as follows: 

Theorem 3.1. Let d G IR 3 oe an arbitrary constant unit vector. There exists two positive 
constants r\ and V, such that, if po — 1 G B^'^, uo G -B 5 and do — d G satisfy 

\\po — 111 - 1 3 + lluoll l + lido — dll - i 3 < 77, 
i/ien system (ll.lj) - (11.2p /ias a unique global strong solution (p, u, d) urai/i (p — 1, u, d — d) 

3 

m OS 2 satisfying 

|| (p - i,u,d - d)!!^ < r (||po - iy,3 + ||u |y + ||d - d||^ H ) . 

3 

The local existence in Q5y can be established by a standard fixed point argument, for 
instance, see [3]. The global existence of fjl . 1[> — HI .2[> will be established by extending a 
local solution with the help of uniform estimates for the local solution when the initial 
data is sufficiently "small". In the rest of this paper, we focus on the uniform estimates 
and uniqueness of the solution to (|l.ip - (|1.2p . 

Remark 3.1. Similar results for other dimensions are still true. As the density is equal to 
one, we recover a global existence result for incompressible liquid crystal flows, which is 
similar to 1231 . 
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4. Uniform Estimates for Linearized Systems with Convection 

In this section, our goal is to obtain uniform estimates of local solutions. For this 
purpose, we consider proper linearized systems with convection that are associated with 
the reformulated system (|3,3|) . First, we investigate the following linearized equations for 
O and d. 

d t n + u • Vn - /xAfi = £, (4.1a) 
d t d + u- Vd- Ad = M, (4.1b) 

where C, A4 and u are given functions. For the system (|4.ip . we have the following 
estimate: 

Proposition 4.1. Denote 

V(t):=J ||u( S )|| B| d S . (4.2) 
Let (f2, d) be a solution of (|4.ip on [0, T). Then the following estimate holds on [0,T): 

-i3+|in(t)|| i+/ d|d(a)||- B 7 +||fi(s)|| Ads 







t 

X ds\ 



- 1 3 



<Ce cl/ Wf||d ||- 1 3 +||n || i + I e- cv ^(\\C\\ i + 
where C is a universal positive constant. 

Proof. To prove this proposition, we first localize (|4.ip into low and high frequencies ac- 
cording to the Littlewood-Paley decomposition. Then each dyadic block can be estimated 
by using energy method. 

Let (fl, d) be a solution of (|4.ip and K > 0. We introduce the following transformation 
of variables [3] : 

n = e- KV V>n, d = e- KV ^d, C = e- KV ^C, M = e~ KV ^M. 



Applying the operator A q to the system (|4.1|) . we deduce that (AgCl, A„d) satisfies 

(4.3) 



5 t A g O + A g (u • Vfi) - fiAA q 9, = A q C - KV'(t)A q n, 



d t Agd + A q (u ■ Vd) - AAgd = A q M - KV'(t)A q d. 

The proof can be carried out in three steps. 

Step 1: Low Frequencies. Suppose q < and define 

/ g 2 = ||A^||2 2 + ||A g d||| 2 . 



Taking the L -scalar product of the first equation of (|4.3p with A q Q, and the second 
equation with A q d, we obtain the following two identities: 

—HA^niUa + (a,(u • vn)\A q n) + m ||aa 9 o||| 2 = (A q c\A g n) - KV'\\A q nf L2 , (4.4) 
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and 

iA||A,d||| 2 + (A g (u • Vd)|A 9 d) + ||AA 9 d||2 2 = (A q M\A q d) - KV'\\A q d\\ 2 L2 . (4.5) 



Adding (|4.4p and (|4.5p together, we obtain 

l_d 
2dt 

where 



/J + H|AA 9 0||| 2 + ||AA ? d||| 2 = X - KV'(t)fl (4.6) 



* = -(A,(u • VQ)\A q n) - (A,(u • Vd)|A 9 d) + (A g £|A 9 Q) + (A q M\A q d). 
The term * can be estimated by using Lemma |2. II (taking si = S2 = § in (|2.ip ) such that 
|*| < C/ 9 (||A 3 £|| L2 + ||A^|| L2 + 2~ia q V'(\\d\\ x + ||f2|| i)) 

; \ (4.7) 

< C/ 9 (||A 9 £|| i2 + ||A ? ^|| L2 + 2-^(11^11^3 + ||n|| B4 : 
In the last inequality, we used the embedding 

-13 1 

Hence, combining (|4.6p and (|4.7p together, we have for q < 
^/ g 2 + ^||AA^||i 2 + ||AA,d||i 2 



2 



(4.8) 



< Cf q (\\A q C\\ L 2 + ||A^|| L2 +2-5^(1^11-1,3 + \\n\\ Bh )j-KV'(t)f ( 

Step 2: High Frequencies. In this step, we assume q > and set 

f q = ||A^||| 2 + ||AA,d||| 2 

We apply the operator A to the second equation of (|4.3|) . multiply by AA g d and integrate 
over M 3 to yield 

— IIAA^II^ + (AA,(u • Vd)|AA,d) + ||A 2 A 9 d||| 2 
= (AA,Ai|AA,d) - KV'\\kA q d\\\ 2 . 



Adding (14. 4p and (|4.9p together, we have 
ld_ 

with 



f q + mIIAA^H^ + ||A 2 A g d||i 2 = y - KV'(t)f 2 q (4.10) 



y = -(A g (u • Vft)|A g fi) - (AA„(u • Vd)|AA g d) + (A 9 £|A g ft) + (AA 9 .M|AA,d). 
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For y, using Lemma 12. 11 one has 



|3>| < ||A 9 n|| L 2||A,£|| £ 2 + \\AA q d\\ L 2\\AA q M\\ L 2+C2~2a q V'\\Q\\ i\\A q n\\ L2 



1 3 



+ C2~^^^-^a q V'\\d\\ - 1 3 ||AA 9 d|| L 2 (4.11) 
< Cf q (\\A q C\\ L 2 + \\KA q M\\ L 2 + 2-ia q V\\\d\\ BU + 



Hence, combining (|4.10p and (|4.1ip together, we have for q > 
1 d 
2dt Jq 



:/ 2 + /i||AA,n||2 2 + ||A 2 A 9 d||2 2 



< Cf q (\\A q C\\ L 2 + ||AA g ^|| L2 + 2-2a,y / (||d||^i,3 + ||fi|| B i)J - KV'{t)f q . 



(4.12) 



Step 3: Damping Effect. We are now going to show that inequalities (|4.8p and 
()4.12p entail a decay for f2 and d. Denote g q = 2% f q for q G Z. The well-known Bernstein's 
inequality implies that 

C g q < 2 2 2\\A q n\\ L2 +2^ M W||A g d|| L2 < l-g q (4.13) 

for some universal positive constant Co- Therefore, we infer from (|4.8p . ()4.12j) that there 
exists a universal positive constant k such that 

I| 5 2 + Kmin{M}2 2^2 < Ca,g,(||£|| B i +\\M\\- hi +ni|d|y,3 +||0|y)) 

-xy'(t)^. 

(4.14) 

Let 5 > be a small parameter and denote x 2 . = 9 q + From (|4.14p . dividing by Xqi we 
obtain 

^ Xg + Kmin{/,,1}2 2 % < Ca q (\\C\\ B i + \\M\\ M + V'(\\d\\ M + 

-Ky / X< ? + 5^' + (5K2 29 . 
Integrating the above inequality over [0, t] and having 5 tend to 0, we obtain, 

g q (t) + Kmin{//, 1}2 2<? / g q (s)ds 
Jo 

<g q (0) + C [\ q ( S )(\\£\\ B i+\\M\\ s i 4 )ds (4.15) 



+ /V'( 8 )(Co ff ( fl )(||d||-i i3 +\\n\\ B t)-Kg q (s)y s . 

J 



dU5]) implies that 

'a9(«)(l|d||ji,: 

< Ca g ( S )||d|| 5 i ; 3 - C ^2^' 5 W|| Ag d|| L2 + Ca ff (a)||n|| D4 - C K2% \\A g Cl\\ L 2. 



Ca q (s)(\\d\\^4 + \\n\\ B i)-Kg q (s) 
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If we choose K such that KCq > C, then we infer from the fact Ylqez a Q — ^ that 

£ {Ca q (s) (\\d\\- hi + ||n|| fl4 ) - Kg q (s)} < 0. (4.16) 

With the inequality (|4.16p in hand, after summation over Z, we deduce from (|4.15p that 
||fi(t)|| fli + l|d(t)||^,3 + Kmm{», 1} |* (||n( S )|| -5 + \\d(s)\\ m ) ds 



(4.17) 



< G{||no|| fl4 + lldoll -,a + j Q (\\c(s)\\ sh + ||^( s )y,a)d s }. 

This finishes the proof. □ 

Next, we turn to consider the linearized system for g and h: 

d t g + u- Vg + Ah = J, (4.18a) 

d t h + u ■ Vh - vAh - Ag = 1C, (4.18b) 

where J , JC and u are given functions. Note that this system is different from (|4.ip since 
there are stronger couplings between g and h. For the system (|4.18|) . we have the following 
estimates (see a similar result in [3] Proposition 2.3]): 

Proposition 4.2. Let (g,h) be a solution of (|4.18p on [0, T). Then the following estimate 
holds on [0, T): 

\w)\\ § i,i + mm B i + f (ii<?wius,j + ii^wiiaj) <*« 

< Ce^){||„|| sM 4- ||M fli + /V-'">(| W | fl , + Ulrf.,)*}, 
where C is a universal positive constant and V is given by (j4.2|) . 

Proof. The proof is due to the argument in [3], and similar to that in Proposition 14.11 For 
the completeness, below we present a proof that is slightly different from [3]. To this end, 
again we first localize (|4.18p in low and high frequencies according to the Littlewood-Paley 
decomposition. 

Let (g, h) be a solution of ()4.18j) and K > 0. Define 

~ Q = e- KV ^g, h = e- KV ®h, J = e~ KV ^J, t = e~ KV ^K. 

Applying the operator A q to the system (|4.18p . we deduce that (A q h, A q p) satisfies 

8 t A q g + A q {u ■ Vg) + AA q h = A q J - KV'(t)A q g, 
d t A q h + A 9 (u • Vh) - vAA q h - AA q g = A q K - KV'{t)A q h. 
Set 

q = log 2 ( - J . (4.20) 
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Step 1: Low Frequencies. Suppose q < q$. As a result, 2 q < 2 qo < ^. Taking 
the L 2 -scalar product of the first equation of (|4.19p with A q g and the second equation of 
()4. 19[) with A q h, we obtain the following two identities: 

ld_ 

2dt' 
and 



lA^Hia + (A,(u • Vg)\A q g) + (AA q h\A q g) = (A q J\A q g) - KV'\\A q g\\ 2 L2 , (4.21) 



~||^|| La + (A,(u • Vh)\A q h) + v\\AA q h\\ 2 L2 - (AA q g\A q h) 

= {A q K\A q h) - KV'\\A q h\\ 2 L2 . 

Next, we derive an identity involving (AA g p| A q h). For this purpose, we apply A to the 
first equation in (|4.19p and take the L 2 scalar product with A q h, then take the scalar 
product of the second equation in (|4.19p with AA q g. Summing up both equalities, we get 



-(AA q p\A q h) + (A,(u • Vh)\AA q p) + (AA,(u • Vp)\A q h) 
- \\AA q pf L2 + ||AA 9 ^|||a + v(A 2 A q h\AA q p) 
= (AA q J\A q h) + (A q IC\AA q p) - 2KV'(AA q p\A q h) 
Let r be a small constant such that < r < |. We define 



(4.23) 



f 2 = \\A q g\\ 2 L2 + \\A q h\\ 2 L2 - ^(AA q g\A q h). 



The Bernstein's inequality yields 



TV 



(AA q g\A q h) 



TV ~ TV 

< —\\AA q ~g\\ L 2\\A q h\\ L 2 < —2<i\\A q g\\ L 2\\A q h\\ L 2 



<2^(\\A q g" 2 



+ i|A^ni 2 ) 



< -(||A^||| 2 + ||A^||| 2 ) 

which implies that 

/ 2 «||A^||| 2 + ||A ? £|| 2 2 . (4.24) 

Here we note that the universal constant due to the Bernstein's inequality is harmless in 
our estimate and thus assumed to be one for simplicity. 

Multiplying (|423j) by and adding it with (|42l) . (|422jl . we obtain that 



1 d 2 



2dt 
with 



r)||AA (? / l || 2 2 +r||AA^|| 2 2 -i,r(A 2 A g / i |AA^) = X 1 - KV' (t) f 2 , (4.25) 



TV 



X l = -(A q (u ■ Vh)\A q h) - (A q (u ■ Vg)\A q g) + — (A,(u • Vh)\AA q g) 



TV 



+ — (AA q (u ■ Vg)\A q h) + {A q K\A q h) + {A q J\A q ~g) 



TV 



{AA q J\A q h) + (A q K\AA q ~g) 
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13 ^ ^ 3 

Since 4>^'^(q) < §, the assumption q < qo implies that 2~ q ^'* 2 ^~^ < 2 qo . As a conse- 
quence, for Xi, using Lemma 12. II and ()4.24p . we have 

< Ca q 2-iV(\\h\\ Bh \\A q h\\ L2 + \\Q\\ B i\\A q g\\ L2 ) 

+ ^a g V(2-*||AA ff e|| £a ||fc|| fl4 +2-^' 4 W- 1 )||^||^,3||A 9 / l || L2 ) 

+ ||A g ie|| L2 ||A^|| L2 + ||A^|| i2 ||A^|| L2 (4.26) 

Ctu 

+ - § -2*(||A g X:|| L2 ||A^|U 2 + ||A gi 7|| i2 ||A^|| i2 ) 

< Cf q (\\A q t\\ L , + ||A 9l 7|| L2 + 2-%a q V'{\\~ e \\ Bl z + \\h\\ Bh \ 

where C is a universal constant that may depend on v. Besides, due to our choice of r, 
we can conclude that 

|i/r(A 2 A,£|AA,e)| < ^IIAA^II^IIAA^H^ < 4||AA g / t ||| 2 + pAA 3 £|| 2 2 . (4.27) 

Then it easily follows from the Bernstein's inequality that 

2 2 «(||A^|| 2 2 + ||A 9 £|| 2 2 ) « (8 - r)||AA^||| 2 + r||AA^||| 2 - ^t(A 2 A^|AA^). 

Hence, combining (|4.25p and (|4.26p together, we can find a positive universal constant k 
such that 

(4.28) 



dt f q + k2 2 *(||A^|| 2 2 + ||A^|| 2 2 



< Cf q (\\\it\\ L 2 + \\A q j\\ L 2 + 2-2a q v'{\\~e\\ B 1,3 + \\h\\ Bk )) - KV'(t)fl 

Step 2: High Frequencies. Suppose q > qo. We apply the operator A to the first 
equation of (|4.19p . multiply by AA q g and integrate over K 3 to yield 

1|||AA^||| 2 + (AA,(u • Ve)|AA 9 e) + (A 2 A,/i|AA,e) 
= (AA 9l 7|AA^)-Ky'||AA^|| 2 2 . 



Set 



/ 2 = ||AA^|| 2 2 + ^||A g ft|| 2 2 - ^(AA^IA^). 



It easily follow from the Cauchy-Schwarz inequality that 

/ 2 «||A^|| 2 2 + ||AA 9 e1| 2 . 
A linear combination of (|4T22|) . ([^23]) and (|4T29|) yields that 

\j t f q + 2 -\\KA q h\\l 2 + i||AA f ^|| 2 2 - 1(AA^|A ? ^) = y x - KV\t)f q (4.30) 
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with 

y l = -^(A,(u • Vh)\A q h) - (AA 9 (u • Vg)\AA q g) + ^(A q lC\A q h) 
+ (AA q J\AA q g) + -(A q (u ■ Vh)\AA q g) + -(AA 9 (u ■ Vg)\A q h) 
- ±(AA q j\A q h) - ^(A q IC\AA q g). 

For 3^i, the assumption q > qo implies that 2~ q ^'^ 2) < 1, then using Lemma |2. 11 we 
can see that 

\yx\ < ^a q 2-iV\\h\\ a l \\*ah\\v + Ca q 2-^ U ^V'\\Q\\ .1 3 ||AA^|| L2 



" ll—g'-lli- " ii^ii^.fii— <?i 

+ A||A 9 £|| i2 ||A^|| L2 + ||AA gi 7|| L2 ||AA^|| L2 

+ i(||AA f/l 7|| L2 ||A^|| i2 + ||A 9 £|| L2 ||AA^|| L2 ) (4.31) 
+ ^n2-*||AA,e|| La ||h|| Bi +2"^' I W- 1 )||^|| M ,3||A^|| L2 ) 
< Cf q (\\A q iC\\ L 2 + \\AA q J\\ Li + 2-ia q V'(\\Q\\ s M + Nl fli ) 
Since q > qo = log 2 (^), the Bernstein's inequality implies that 

^\{AA q ~g\A q h)\ < 2-»l||AA 9 e|| £a ||AA,^|| La < ^(||AA^||| 2 + ||AA^||| 2 ). 
As a result, 

||AA ? £||| 2 + IIA^Hl. < ||AA^||| 2 +2- 2 ^||AA 9 / i ||i 2 

< ||AA^||| 2 + ^-||AA^||| 2 (432) 

< C 0||AA^||| 2 + i||AA^||| 2 - ^(AA q g\A q h) 

Hence, combining (|4.30|) . (|4.3ip and (|4.32p together, there exists a positive constant k 
such that 

(4.33) 

Step 3: Damping Effect. We now show that inequalities (|4.28p and (|4.33p entail 
a decay for h and g. Denote g q = 2% f q for 9 £ Z. It follows from (|4.28p . ()4.33p . and 
Bernstein's inequality that 

- — (a 2 ) + k2^ 2 '°^ a 2 , < Ca.a J WKW 1 + 11.711 -1 3 +V'(\\d\\-i 3 +\\h\\ i)) 

(4.34) 

iq- 



di f 2 q +K(\\A q h\\l 2 + \\AA q g\\ 2 L2 j 
< Cf q (\\A q K\\v + \\AA q J\\ L2 +2-ia q V'(\\g\\ B is + \\h\\ 1 )) - KV'(t)fl 



KV'(t)g 2 q 
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Let 5 > be a small parameter (which will tend to 0) and denote x q = 9 q + <5 2 - From 
(|4,34p . dividing by x<j, we obtain 

± Xq + K 2^ Xq < c^m^ + \\j\\-^ + v(m\\ 6 m + ii^iy: 

Integrating the above inequality over [0, t] and having 5 tend to 0, we obtain, 



t r-t 



g q (t) + K2^°M J % ( 8 )da<^(0) + C7 J a g ( S )(||^|| B i +11^11^1,3)^ 

+ jf ^(Ca,( a )(||e|| 5jiS + \\h\\ B i)- Kg q (s))ds. 



(4.35) 



Bernstein's inequality implies 

C g q < 2§ ||A^|| L2 +2^' 4 W||A g/ 5|| L2 < 

for some universal positive constant Co, and hence 

Ca q (8)(\\g\\ M + fh\\ B i) - Kg q (s) < Ca q (s)\\h\\ fl , - C ^2§ || A^|| L2 

+ Ca q (s)\\e\\~ h s -C K2^ h4 ^\\A q g\\ L ,. 

If we choose .fTCo > C, we have 

2 {c« 9 (s) (||£iy ,3 + ||^|| B i) - #<fo( S )} < 0. (4.36) 

With the inequality (|4.36p in hand, after summation over Z, we deduce from ()4.35p that 

wkm^ + \\m\\ S M + « f (n%)iUu + n^)iy,i) ^ 

° t (4-37) 

< c{iiM fli + ikoiy,3 + y o + njwii fi4lS )d-}. 

Step 4: Smoothing Effect of h. Based on the damping effect for g, we can now 
further get the smoothing effect of h by considering (|4.18p with being seen as a source 
term. Indeed, thanks to ()4.37p . it suffices to state the proof for high frequencies only. We 
therefore assume that q > qo- 

Define I q = 22\\A q h\\ L 2. Then, from the energy estimates for the system 

d t A q h - vAA q h = -A q (u ■ Vh) + AA q g + A q IC - KV'(t)A q h, 

we have 
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for a universal positive constant k. Using = 1^ + <5 2 , integrating over [0,t] and then 
taking the limit as 5 — > 0, we deduce 

I q (t) + v2 2 « [ I q (s)ds < I q (0) + f 2 q i\\A q lC{s)\\ L 2ds+ [ 2%\\A q Q(s)\\ L 2ds 

Jo Jo Jo (438) 

+ cj o V'(s)a q (s)fh(s)\\ Bh ds. 

We therefore get 

£ 2i\\A q h(t)\\ L2 +u f 2%\\A q h(s)\\ L2 ds 
<?><?o i>qo 

<IIM i+ l|/C( S )|| ida+ fE 2f IIA^WH^cfo + CY(t) sup ||/i|| i. 
Using (j4.37|) . we eventually conclude that 

"fE 2^||A^( 8 )|| ia cfa < C(l + U(i))(||^o||^i,3 + 

+ r(^wu+ii^wiUi.i)*)- 

u 

Combining the last inequality with (|4.37p . we finish the proof of Proposition 14.21 □ 



5. Global Existence for Initial Data Near Equilibrium 
In this section, we are going to show that if the initial data 

llpo — 111 - i 3 + II Un II 1 + lldn — d|| - i 3 < n 
for some sufficiently small rj, there exists a positive constant V such that 

H( P -i,u,d-d)y <r v . 

This uniform estimate will enable us to extend the local solution (p, u, d) obtained within 
an iterative scheme as in [3] to be a global one. To this end, we use a contradiction 
argument. Define 

r = sup{TG [0,oo) : ||(p - 1, u, d - d)|| s < Tr/}, 

with r to be determined later. Suppose that Tq < oo. We apply the linear estimates in 
Proposition 14. II and Proposition 14.21 to the solution of reformulated system (|3.3p such that 
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for all t E [0,To], the following estimate holds: 

rT 



< Ce cv (\\d - d|| -i a + ||Oo|| x + ||£|| rl + 



I 3 



(5.1) 



and 

< ^ (||po - 1|| -,3 + ||M fli + m\ L}roiBh + 11^11^(^,1)) > 

where 

fTo 

V = I Hull 5 ds. 



(5.2) 



'B5 

We note that Vd = V(d — d) because d G S 2 is a constant vector. As a result, the 
function A/" (see (|3.1f> ) can be rewritten as 

M = M := -u • Vu - ^— L4u - -div ( V(d - d) V(d - d) - -|V(d - d)| 2 / 

p p V 2 

Therefore, the functions C and Ai in (|5.ip are given by 

C = A -1 cur]A? + u- Vfi, 

a^ = |v(d - d)| 2 (d - d) + |v(d - d)| 2 d, 

while for the functions tC and J in (|5,2p . we have 

/C = A~ 1 divAA + u • V/i, 
J = -{p- l)divu. 

In what follows, we derive estimates for the nonlinear terms C, A4, tC and Indeed, 
for the term u • Vf2, we infer from Proposition 12.21 that 

/•To 

<c||u|| i ii^H . 5 <crV- 

Similarly, we have 

rT 



||A" Wu • Vu|| fl4) < C||u • Vu|| j } < ^ ° ||u|| flJ ||Vu|| s| 



3 dt 



< C\\u\\ , i llull . 5, < CT 2 r] 2 . 



(5.4) 
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By the embedding <—> L°°, we consider F(g) = j^- with q = p — 1, then by [H Lemma 
1.6], we have 

\\m\\ B i <c*{\\q\\l~)\\q\\ b z, 

where ||£>||l°° < 1 + ||p||l°° < 1 + C||p|| a . Then we can apply Proposition 12.21 to obtain 

B ^ 

that 



A _1 curl Qdiv ^V(d - d) V(d - d) - i|V(d - d)| 2 /^ 



<C 1 + 



p-1 



V(d - d) V(d - d) - i|V(d - d)| 2 / 



<C l + ||p-l|| 



L»(B*) 



||v(d - d) v(d - d) - -i v(d - d)| 2 /||^ {b1) 



<c i + lip -in -i3 N 



To 



\\V(d-d)\\lsdt 



\ 2 - 



To 



<C(1 + Tri) I ||V(d-d)||^3dt 



< C(l + rr?)||d-d| 



<C(l + rn)||d-d|| -islld-dll -57 

< c(i + rr ? )r 2 r ? 2 . 



In the last line of above inequality we have used the interpolation 



(5.5) 



8 « <C|I/II 



The term ^-y-Au can be dealt with in a similar way such that 



A x curl 



p-1 



Aw 



< 



P-1 



< 



P 

P-1 



A 



u 



P 



\\Au\\. 



<C\\p-l 
<C\\p-l 
< CT 2 ri 2 . 



- 1 3 , llull , 5, 



(5.6) 



Combining (I5.3p - (l5.6p together, we obtain 



\\c\\ L ^ {Bh) <cr 2 v 2 + cr 3 v 3 . 



(5.7) 
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Similarly, we can derive the estimate for /C: 

||ic|| i < cry + cry. (5.8) 

Next, we turn to the estimate for J\4. For |V(d — d)| 2 (d — d), indeed, we have, as in 
|||V(d-d)| 2 (d-d)||^ (5j(J) 

^GIIIvCd-ajHi^iKd-a)^^,,, 

/•To 

<CT^ / ||V(d-d)|| 2 3d* 
Jo B1 - 

<Cr ?? / ||V(d-d)|| 2 , ,dt 

Jo B2 2 

<c^ii(d-d)n^ ( . H) 

^^iKd-^n^^iKd- ^\ Lk{§U) 

< CT 3 t] 3 . 

For |V(d — d)| 2 d, we have, by the definition of Besov's spaces 

lll v ( d - d )l% o (BH)^lll v ( d - d )l 2 ^ o( ^) 

(■T 



<c j\\v{a-a)\\ L ~\\v(a-d)\\^ u dt 
<c /" T °||v(d-d)|| 2 i 3 dt 

Jo B2 2 

< C||d-d|| 2 a 5 X 

< C||d-d]| -i s Jld-dll . -5 7 N 



< cry. 

Finally, for we have 



||(p — l)divu|| , -i3 < Clip— 111 -i3 l|divu|| . 3 N 

llU ' "Ll, o (B2'2) - l|p M L°° (B3'I)" M Ll, (S3) 



(5.9) 



(5.10) 



(5.11) 



< cry. 

Substituting ([577 ]) — ([57TT]1 back to d57T]) and ([572]), we obtain 

||( P -i,u,d-d))|| a < c ie ^ r " („ + ry + ry ) . (512) 



SB- 
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We choose T = 4Ci, and then r\ > satisfying 

e^<2, r%<i rV<i (5.13) 
Hence, it follows from (|5.12p and the above choices of T and r] that 

||(p-l,u,d-d)|| 3 <T V , 

which is a contradiction with the definition of To. As a consequence, we can conclude that 
Tq = oo. The proof of global existence is thus proved. 



6. Uniqueness 

3 

In this section, we will address the uniqueness of the solution in . For this purpose, 

3 

suppose that (pi, Uj, dj)j = i ; 2 in 55 2 solve (jl.ip with the same initial data. 
Define 

Sp = p2 — pi, Su = U2 — ui, 5h = A _1 div5u, <5f2 = A _1 curl5u, 5d = d2 — di. 

3 

then (5p, 5u, 5d) £ *8y for all T > 0. On the other hand, since (pi, Uj, dj)j = i 2 ar e solutions 
to (jl.ip with the same initial data, (5p, 5u, 5d) solves 



d t 5p + u 2 • V5p + A5h = 5 J, 

d t 5h + u 2 • V5h - vA5h - A5p = SK., 

d t sn + u 2 • V5n - A5n = sc, 

d t Sd + u 2 • VSd - A5d = SM, 
Su = -A^VSh + A-^urWO, 



(6.1) 



with 



SK 



5 J = —Su ■ Vpi — (5/)divu2 — (p± — l)div5u, 
u 2 • VSh + A _1 divf - u 2 • V<5u - Su ■ Vui - ( ) ^lu 2 + — ASu 

\ \Pl P2j Pi 

- — div(V(d 2 - d) VSd + VSd V(d x - d) 
Pi v 

- -/(V(d 2 - d) : VSd + VSd : V(di - d 

— -—) div(v(d 2 - d) v(d 2 - d) - i|v(d 2 - d)| 2 / 

Pi P2j v 2 
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5C = u 2 • V5Q + A _x curl( - u 2 • V<5u - Su ■ Vui - I ) Au 2 + —ASu 

v \pi P2J pi 

- — div(V(d 2 - d) V5d + V5d V(di - d) 
Pi v 

- ^/(V(d 2 - d) : V5d + V5d : V(di - d) 

1 _ L\ div(v(d 2 - d) © v(d 2 - d) - i|v(d 2 - d)| 2 / 

Pl P2j v 2 



1 



and 



= |V(di - d)\ 2 5d + (V£d : V(di - d) + V(d 2 - d) : V5d)d 2 , 



where we used the notation A: B = Ylij=i ^-ij^ij- 

Applying Proposition 14. II and Proposition 14.21 to the system (|6.ip . we get 



\\(Sp,Su,Sd)\\ f 

< exp (c\\M^ Bi) ) (ll^ll^Lij + ll^ll^^-l) + II^ ( b-4) (6-2) 

+ 11^(5-1.4, 

Since 5u, <5d) G ®|> — 1) € L} oc (B%), and hence # - 1 G C(Si) n L°°(sf). 
This entails — 1 G C([0, oo) x R 3 ). On the other hand, if 77 is sufficiently small, we have 

\pl(t,x) - 1| < - for all £ > and x G I 3 . 
Continuity in time for p 2 — 1 thus yields the existence of a time T > such that 

- 1||l°° < x for * = 1, 2 and * € [0, T]. 
Repeating the argument in Section 5, we easily infer that 

\\$J\\ 1 - ,-1 sjtfull , 3 N + ||divu 2 || 3 \\5p\\ - ii, 

~ Ii ua iii4(fli)ii v * u iixi.(fl-i) + ii 5u II^(^)II Vui IIl| ( b^) 



+ ( 1 + llPi - lll^^J, + II/* - lll^^ijj II V'u 2 ||^ (Bi) ||5p|| LS?(si) 



+ ( 1 + llPi - l|l L . (Bi) + II*. - lH^j |||V(d 2 - d)|i^ iBi) \\8p\\ h) 



+ (1+11/91 -111 , 3j||M|| , -13 (||di-d|| , -3 5 + ||d 2 - dll , -3 s N 



and 
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\SM\\ , - 11 <||<5d|| - i i ||V(di -d)|| 2 ,3 5 



l + ||d 2 -d|| , s J ||V(d 2 -d)|| . , a, 



+ ||V(di -d)|| . , 3 J||V(Jd|| - i i . 
Substituting those estimates back into (|6.2p . we eventually get 



||(^,<Ju,(5d)|| i <Z(T)||(<^u,Sd)|| 

with 



Z(T) = e X p[C\\u 2 \\ LUBh 



I Pi — 111 ,-i3+ 1 1 divuo 1 1 , , a. + 1 1 u-2 1 1 „ , a, 



+ ( 1+ ^- 1 "W(-)) l|v(dl - a)ll ^ ) 



+ i + IIpi - ilLoc^a, + IIP2 - i|| roo , 3 J|||v(d 2 - d 



|2| 



3 



+ I 1 + lIPl — 111 , 3 , ) ( ||dl — dll „ - 3 5 % + ||do — dll „ -35, 



+ I 1 + "" d2 - d h ¥i Bh) U" V ^ - d A^h + l|v(dl " d)ll ^ ( ^) 

We notice that limsup r ^ + Z(T) < C\\pi — M\ Loo \- This is because all other terms 
involve an integral in time in L 1 or L 2 sense so that as T goes to zero, all those integrals 
will converge to zero. Thus, if 7] > is sufficiently small, we get 

||(*MMd)||^ (fl4) = o 

for certain T > small enough. Thus, we have shown the uniqueness on a small time 
interval [0, T] such that (/?i,ui,di) = (p2,U2,d2). 

Then we can argue as in [3] for the compressible Navier— Stokes equations. Let T max < 
+oo be the largest time such that the two solutions coincide on [0, T max ]. Taking T max as 
the initial time, we denote 

def 

{pi(t),Ui(t),di(t)) = {pi(t -T max ),Ui(t -T max ),di(t -T max )). 

Repeating the above arguments and using the fact that ||pi(0) — < j, we can prove 

that 

(pi(t),ui(t),d7(t)) = (^(t),u5(t),$(t)) 

on a sufficiently small interval [0, t] with i > 0. This contradicts the assumption that T max 
is the largest time such that the two solutions coincide. Thus, T max = +oo which means 
that the uniqueness result holds in ]R + . 
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